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ABSTRACT 



Aims. To estimate the feasibility of dense-matter phase transition, we studied the evolution of the central density as 
well as the baryon chemical potential of accreting neutron stars. We compared the thin-disk accretion with and without 
the magnetic field torque with the spin-down scenario for a selection of recent equations of state. 

Methods. We compared the prevalent (in the recycled-pulsar context) Keplerian thin-disk model, in which the matter 
is accreted from the marginally-stable circular orbit, with the recent magnetic-torque model that takes into account 
the influence of stellar magnetic field on the effective inner boundary of the disk. Calculations were performed using 
a multi-domain spectral methods code in the framework of General Relativity. We considered three equations of state 
consistent with the recently measured mass of PSR J1614-2230, 1.97di0.04 M0 (one of them softened by the appearance 
of hyperons) . 

Results. If there is no magnetic torque and efficient angular momentum transfer from the disk to the star, substantial 
central compression is limited to the region of initial stellar masses close to the maximum mass. Outside the maximum 
mass vicinity, accretion-induced central compression is significant only if the angular momentum transfer is inefficient. 
Accounting for the magnetic field effectively decreases the efficiency of angular momentum transfer and implies a 
significant central compression. 

Conclusions. An efficient angular momentum transfer from a thin disk onto a non-magnetized neutron star does not 
provide a good mechanism for the central compression and possible phase transition. Substantial central compression 
is possible for a broad range of masses of slowly- rotating initial configurations for magnetized neutron stars. Accretion- 
induced central compression is particularly strong for stiff equation of state with a high-density softening. 

Key words, dense matter - equation of state - stars: neutron - stars: rotation 



1. Introduction 

Neutron stars (NSs) are detected in binary systems by X- 
and 7-ray observatories because of the immense amounts 
of radiation released in the process of accretion of matter 
falling from a companion star onto the NS surface. Because 
the stellar mass increases during accretion, it is usually as- 
sumed that the density in the center of the star increases 
as well, and that a sufficient increase of mass may result 
, in a phase transition (the details of which are not yet fully 
. known). It is expected that a phase transition in the NS 
core might cause observable astrophysical phenomena, such 
, as star-quakes, spin clustering, gravitational wave emission 
■ and therefore serve, it is hoped, as a testbed for the quan- 
' tum c hromodynamics (QCD) phase d iagram (ICheng fc Pail 
" 19981: iGlendenning fc Webed l200ll: iBlaschke et all I2OOII: 



Zdunik et al1l2QQ5l : iBlaschke et al.ll2Q"Q8[ IXu fc Liang! l2QQ9l 

and references therein). 

We aim to show with three recent equations of state 
(EOSs) of NS cores, that the central baryon number den- 
sity, ric, and the baryon chemical potential, /ic, may change 
in a non-trivial way while the star gains in mass and si- 
multaneously increases its rotation rate during accretion. 
In particular, we will study whether, and if so, to what 
extent, the evolution of ric and fic depends on the EOS 
of dense matter. Our considered EOSs are consistent with 
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the rece ntly-measured niass o f PSR J1614-2230, 1.97 ± 
0.04 M0 (|Demorest et al.L [20101 ) . Complementary questions 
that may be asked are (a) how does the actual efficiency of 
the angular-momentum transfer and the stellar magnetic 
field coupling to the disk affect the final rapidly-rotating 
NS configuration, and (b) how does the accretion-induced 
evolution depend on the mass of the initially slowly rotating 
configuration. Thus far, the formation of millisecond pul- 
sars was most ofte n modeled following the classical paper of 
ICook et al.l (|l994l ) who assumed the stellar magnetic field 
is sufficiently weak that it does not influence the accretion 
flow significantly; consequently, one assumed B = and the 
accretion from the marginally-stabl e circular orbit (of ra - 
dius Vms)' No t able e x ceptions were iPossenti et al.l (119991), 
iBurderi et aD (I1999D. IBlaschke et al.l (|200l[ ). IColpi et al.l 
(120011) . and'Blas chke et all (|2008f ), who used simple models 
of magnetic torq ue acting in the NS-accretion-disk system 
(iGhosh fc Lama I1979L see also IGhosh fc Lambl[T99ll ). 

Here we compare the results obtained using the B = 
approach with those employing the recent e f fectiv e 
magnetic-torque model of iKluzniak fc RappaportI (|2007[ ). 
used recently by iBeiger et al.l (|2011 ) to study the forma- 
tion of a millisecond pulsar PSR 1903+0327. The details 
of magnetic-fleld decay during the accretion-driven evolu- 
tion phase of an NS' life are still far from being under- 
stood in spite of theoretica l and observationa l efforts in 
recent years (see Sect. 4 of iBeiger et al.l I2OIII for a brief 
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summary). Hence, we employ the simplest possible descrip- 
tion of the magnetic-field decay that was shown to be con- 
sistent with at least a subclass of observed accreting NSs 
(|Shibazaki et al.L Il989f ) - this approach allows us to study 
the qualitative features of the influence of disk accretion on 
the behavior of the NS central parameters. 

The text is organized in the following way. Sect.[2]gives a 
short description of the methods used to obtain the accret- 
ing NS tracks. The results are presented in Sect. [3l Sect. |4] 
contains our summary and conclusions. 



2. Methods and models 

We simulated the spin evolution of NSs caused by the 
thin-disk accretion without a magnetic field and com- 
pa red it with the results f or the magnetic-torque model 
of iKluzniak fc RappaportI (|2QQ7[ ). Because the original 
Newtonian magnetic-torque model was insufficient in ex- 
treme cases of rapidly-rotating, massive and compact NSs, 
we ad ditionally employed the modification by Beiger et al.l 
(|2Qll[ ), which takes into account the relativistic effects 
caused by the existence of the marginally-stable circular 
orbit (recalled briefly in the appen dix). We considere d 
three recent EOSs - APR EOS by lAkmal et al.l (|1998[ ). 
DH EOS by iDouchin fc Haensell (|200lD and BM EOS 
from a set of non-l i near relativistic mean field models of 
iBednarek fc Mahkal (|2QQ9li n. The first two EOSs assume nu- 
cleon NS cores, while the BM EOS has a characteristic high- 
density softening associated with appearance of hyperons. 
Constructing an EOS yielding M^ax > 1-97 in spite 
of the hyperon softening is not an easy task, therefore it 
deserves some additional explanations. The crucial feature 
of our BM EOS is a repulsive contribution to the pressure 
coming from the quar t ic terms in the vector meson fields 
in the Lagrangian (.Bednarek fc Mahkal . 120091 ) . We empha- 
size, that the BM model reproduces semi-empirical nuclear- 
matter data as well as the semi-empirical estimates of the 
potential wells of hyperons in the nuclear matter, com- 
ing from hypernuclear data, involving hypernuclei and E~- 
atoms, as well as the strength of the A — A a ttraction in the 
AA- hypernuclei (|Bednarek fc Mahkal . l2009t ). Neutron-star 
configurations were obtained with the numeri cal library 
LORENE0, using the implementation of the Bona zzola et al.l 
(|1993[ ) formalism for axi-symmetric and rigidly-rotating 
stars (rot star code). 



2.1. Magnetic torque neglected 

Assuming that the magnetic field is not disturbing the 
Keplerian thin accretion disk structure, one usually con- 
siders the spin-up scenario in which the angular momen- 
tum is transferred from the disk to the NS via accretion 
from the mar ginally -stable circular orbit (|Cook et al.Lll994l : 
IZdunik et ani2002[ ). In this process, an infall of a particle of 



according to 



baryon mass dMb and specific orbital angular momentum 
/ms leads to a new quasi-stationary stellar configuration of 
baryon mass Mb + dMb and angular momentum J + dJ, 



^ We use a specific model corresponding to the parameter 
Av = 0.016. It fits semi-empirical hypernuclear and nuclear data 
and yields for an NS rotating at 317 Hz (spin frequency of PSR 
J1614-2230) Mn,ax = 2.01 Mq. 

^ |http : //www . lorene . obspm . f r | 



dJ 

dMb 



(1) 



where the parameter xi (< 1) quantifies our lack of knowl- 
edge of the fraction of angular momentum that is trans- 
ferred to the star by an infalling particle; recent numerical 
simulatio ns suggest that the v alue of xi is close to u nity 
(see, e.g.. lBeckwith et al.|[2QQ8l and lShafee et al.ll2008D . 

2.2. Magnetic torque included 

If the magnetic field does affect the accretion flow, the 
B = model ceases to be correct. To calculate the in- 
crease of the total stellar angular momentum J in this 
case, one must modify the Eq. ([T]) to incorporate the fact 
that the disk now terminates at some ro > r^^, depend- 
in g on the value of magne t ic fiel d. We use the prescription 
of iKluzniak fc RappaportI (|2007[ ) , which gives the following 
evolution equation: 



dJ 

dMb 



9r3Mb 



3-2 



^0 



3/2' 



(2) 



where /i = BR^ is the dipole magnetic moment of an NS 
(assumed to be parallel to J), rcor is the corotation radius 
and Mb denotes the mean accretion rate; for a detailed de- 
scription of the parameters and modifications in obtaining 
ro, related to the relativi stic marginally-stab le circular or- 
bit, see the appendix and lBeiger et al.l (|201ll ). 



1600 
1400 
1200 
_ 1000 
— 800 
600 
400 
200 




mass-shedding 
instability 
spin-up {xi=l) 
spin-down , . ♦ 




Fig. 1. (Color online) Evolution of the central baryon num- 
ber density ric in the disk-accretion spin-up as well as the 
spin- down scenario, shown for the spin frequency /-central 
baryon density ric plane (DH EOS, B = 0). 

A complete evolution model of a magnetized NS ac- 
creting from a disk needs a formula for the accretion- 
induced surface magnetic field decay. Because the the- 
oretical state-of-the-art of this subject is quite com- 
plicated and far from being complete (see Sect. 4 of 
'Bei ger et al.l l201ll for a brief summary) , we used the sim- 



bvlShibazaki et al.l (Il989|: see also 


Taam 


fc van den Heuvell 


1 19861: Ivan den Heuvel fc Bitzaraki 


19951) 


and assumed that 



B decreases as a function of accreted mass AM only: 
B = Bi/{1 + AM/tub), where tub = 10~^ M©. We also 
tested other available forms of decay law (exponential-like 
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formula used recently bv lOslowski et al.ll2Qlll in the context 
of population-synthesis studies, and the quadratic modifi- 
cation of Shibazaki et al. law) and conclude that this choice 
has no qualitative influence on the results. 

3. Results 

To estimate the feasibility of dense-matter phase transi- 
tions, we examined the compression of matter in the cores 
of accreting NSs by measuring the increase of the central 
baryon number density ric as well as the central baryon 
chemical potential /ic = (Pc + ec)/nc {pc and Cc are the pres- 
sure and mass-energy density, respectively). Furthermore, 
we compared the results for spin-up-induced central com- 
pression with the isolated NS spin-down results. 

3.1. Magnetic torque neglected 

For B = the features of accretion-driven evolution of Uc 
and /ic are qualitatively similar for all investigated EOSs. 
We will thus illustrate these features using the DH EOS 




[fm-^] 



Fig. 2. (Color online) Accretion tracks of the gravita- 
tional mass-central baryon number density, M — nc, plane. 
Magnetic field tracks were obtained for Mb = 10~^ Mq/jt. 
The shaded region contains configurations for which the 
stellar equatorial radius i?eq is larger than r^s (DH EOS). 

and begin with the comparison of spin-up (with the most 
efficient angular momentum transfer xi = 1) and spin- 
down-induced increase of Uc pictured in Fig. [T] in the spin- 
frequency-central density plot; the configurations are al- 
lowed to span the whole frequency range from the initially 
non-rotating stars up to the mass-shedding limit (or the axi- 
symmetric perturbation instability limit). Evidently, the 
spin-down allows for a higher ric increase than the xi = 1 
accretion spin-up. To investigate this feature. Fig. [2] shows 
the behavior of the gravitational mass M as a function of 
Tie. For an astrophysically- motivated range of NS initial 
masses. Mi = 1.2 — 1.6 M©, the xi = 1 accretion never 
leads to a substantial central compression; moreover, the 
maximal central compression is attained for moderate spin 
periods, followed by a decrease of Uc. On the contrary, low- 
ering the efficiency of angular momentum transfer to e.g. 
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Fig. 3. (Color online) Maximum central density nc,f attain- 
able for a given evolutionary track compared to the initial 
nc,i, both plotted as functions of the gravitational mass of 
the initial non-rotating configuration, M[. Upper panel: 
spin-up scenario. The reference line nc,i is the static se- 
quence. Lower panel: spin-down scenario. The reference 
nc,i Hue is the mass-shedding limit, to be compared with 
nc,f of static configurations (DH EOS). 

xi = 0.5 allows for a substantial Uc increase. A detailed 
summary and comparison between various types of tracks 
is included in Fig. [3l where the initial and final Uc is plotted 
as a function of the initial stellar mass Mi for both spin- 
up (upper panel) and spin-down (lower panel). The nearly 
horizontal line for xi = 0.5 in the upper panel denotes the 
instability limit. Based on these data we conclude that to 
achieve a significant Uc increase in an astrophysically-sound 
evolutionary process one should consider 

• an initial mass of the configuration quite close to the 
value of the static-star maximum mass in case of xi = 1, 

• or a value of xi significantly lower than 1 (dotted lines 
in Fig. [2] and [3]), 

• a more efficient mechanism, such as the spin- down of 
an isolated NS - constant Mb track leads to a higher 
increase of Uc than thin-disk accretion, as shown in Fig. 
[T] and the lower panel of Fig. [3l 

Results for the chemical potential /ic for the DH EOS are 
plotted in Fig. IH results for APR and BM EOSs are pre- 
sented in Figs. [Hand [6l 

3.2. Magnetic torque included 

An inclusion of the magnetic torque acting in the NS-disk 
system changes the results significantly, as Fig. [2] shows. 
As an example, we used the "canonical case" of the initial 
magnetic field Bi = 10^^ G and the average accretion rate 
Mb = 10~^ Mq/jt (dashed lines; dots denote the decaying 
magnetic field for the initial mass Mi = 1.4 M©). Overall, 
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the magnetic-torque results resemble those for a reduced 
efficiency of the angular momentum transfer, xi = 0.5. 
Including the magnetic field substantially increases Uc and 
/ic for an astrophysically interesting range of initial masses, 
Mi = 1.2 - 1.6 Mo (Figs. [3] and 11 for the DH EOS, and 
Figs. El and [6] for APR and BM EOSs). The nearly hor- 
izontal line segments in the upper panels of Figs. [3] and 
[5] correspond to the axisymmetric instability limit - it is 
reached by configurations that were not spun up to the 
mass-shedding limit; lowering the mass accretion rate re- 
sults in an extended instability limit line toward smaller 
initial masses. In addition to the maximum attainable com- 
pression, the results for a predefined amount of accreted 
mass (AM = 0.1 M©, AM = 0.25 Mq) are plotted for 
comparison in Fig. [3l 

For the DH EOS and M = 1.4 M©, accretion of about 
0.25 M0 implies (nc,f — ^c,i)/^c,i ^ 0.2. This should be 
compared with very small (of a few percent) compression 
for 5 = and xi = 1 case. Including the magnetic torque is 
therefore crucial for accretion-induced compression in the 
core of an NS with initial mass ^ 1.4 Mq. Note also a 
fairly strong effect of the EOS stiffness below approximately 
three nuclear densities, as well as its high-density behav- 
ior. As seen for example in Fig. [5l for our fixed magnetic 
dissipation evolution model, a stiff BM EOS with a high- 
density hyperon softening predicts a fractional increase of 
Tie by about 120% at Mi = 1.4 Mq , while /ic increases 
by ~ 500 MeV. The overall picture obtained for magne- 
tized stars, with Mi = 1.2 — 1.6 M©, is therefore essentially 
different from the B = one. 
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Fig. 4. (Color online) Gravitational mass of the initially 
non-rotating configuration. Mi as a function of the maxi- 
mum attainable central baryon chemical potential /ic,f for a 
given evolutionary track compared to the initial /ic,i- Lines 
are denoted analogously to Fig.[3l Left panel: spin-up sce- 
nario. Right panel: spin-down scenario (DH EOS). 



4. Conclusions 

We have studied the influence of two thin-disk accretion 
models - with and without the involvement of magnetic 
field - on the evolution of the central density and chemi- 
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Fig. 5. (Color online) Increase of the central baryon num- 
ber density tIq as a function of the initial gravitational mass 
Mi for the APR and BM EOSs. Notations as in Fig. [3l 
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Fig. 6. (Color online) initial gravitational mass Mi as a 
function of intial and maximal central chemical potentials, 
/ic for the APR and BM EOSs. Notations as in Fig. H 

cal potential in accreting NSs. For B = accretion from 
the marginally-stable circular orbit leads to a negligible 
central compression (a few per cent) for astrophysically- 
relevant initial configurations (slowly-rotating stars with 
Mi = 1.2 — 1.6 M0) and maximally-efficient angular mo- 
mentum transfer {xi = 1 in Eq. (TJ a value currently favored 
by numerical simulations). Consequently, this mechanism 
does not seem viable to trigger a dense-matter phase transi- 
tion. Substantial central compression could be obtained for 
xi smaller than unity, or in a different evolutionary process 
(e.g., spin-down). Including the magnetic torque changes 
the outcome profoundly, but for a typical initial magnetic 
field, Bi 10^^ G and accretion rate Mb = 10~^ Mq one is 
able to produce massive, weakly-magnetized {B 10^ G) 
millisecond pulsars with a sizable (even as large as 100%) 
central compression, thus probing much higher densities, 
which in turn may trigger a phase transition. Even such a 
low final value of B still influences the accretion process 
and should be accounted for; only lowering it even more, to 



Bejger et al.: Evolution of central parameters in accreting NSs 



5 



~ 10^ G, essentially allows for the recovery of the B = 
results. 

Note also a quite impressive dependence of the NS spin- 
up evolution on the very presence of the magnetic-field. The 
magnetic torque substantially decreases the spin-up effi- 
ciency, as the comparison with the xi = 1 and xi = 0.5 
results show - in other words, spinning- up with magnetic 
field requires more accreted mass to reach a desired spin fre- 
quency. A considerable subset of " magnetic-torque" tracks 
ends at the axisymmetric instability limit (i.e., nearly hori- 
zontal lines in Figs. [3] and [5]), while the vast majority of the 
xi = 1 tracks can reach the mass-shedding limit. This result 
correlates with the EOS stiffness and is most pronounced 
for the BM EOS, softened by the existence of hyperons; 
potentially, it may help in the understanding of a puzzling 
non-detection of submillisecond pulsars, as well as assist in 
the studies of formation of stellar-mass BHs and their mass 
function. 

We have restricted ourselves to the EOSs of hadronic 
matter. A phase transition to quark matter softens the 
high-density EOS, so that reaching M^ax > 1-97 Mq 
requires some tuning of both hadronic and the quark- 
matter model. An approach based on an effective 
model of the QCD of qua rk matter used recently by 
iBonanno fc SedrakianI (|2Qllf ) indicates that the vector re- 
pulsion in quark matter should be sufficiently strong to 
reach M^ax > 1-97 Mq (we recall that vector- meson re- 
pulsion in hadronic matter is also crucial for our BM EOS 
model). In any case, the maximum mass of NSs with quark 
cores (so-called hybrid stars) turns out to be very close to 
that reached at the central density equal to the deconfine- 
ment density ([Bonanno SedrakianI . 1201 ID . 

Our results were obtained using three specific EOSs of 
dense matter, an effective magnetic torque model in the 
pulsar-accretion disk system, and a simplistic description 
of accretion-induced magnetic-field dissipation. We believe, 
however, that in a qualitative sense these results posses a 
general validity. 
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Fig. A.l. (Color online) Difference (in km) between the 

real value of rms and the value rml*^^"^ (to be compared with 
average radii of the problem: i^ea, ^ ms ^ 10 kp) obtaine d 
by means of an approximation of Bej ger et al.l (|20101[2QTi[ ), 
plotted as iso-contours on the M- f plane (DH EOS). 

near rms 

_ l-Q;/f3/2 /f-2 v'^ 3a \ 

(l-i;2/c2)3/2yrTrT7^ (v^"^^^r3/2-ay' ' ^ '^^ 

with a = Jc/{V^GM'^), f = vq/vs and Vs = 2GM/c^ denoting the 
Schwarzschild radius. Marginally stable circul ar or bit radius rms is 
then a solution of equation /ms(^ms) = 0. Fig. lA.ll shows, using the 
example of the DH EOS (other EOSs we tested yield similar results) 
the comparison of the rms calculated exactly from integrals of motion 
with the approximate value rmf \ obtained by means of Eq. (|A.2|) . 
For almost all configurations the difference is less than 0.6 km (for 

rms < R we adopt rms = R)-, and the regions where |rms — ^ms*^^ '^l > 
0.6 km are small and well-confined near high-mass and mass-shedding, 
submillisecond rotation. 
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Appendix A: Calculation of the inner-boundary 
radius of the accretion disk 

To simphfy the calculation of the specific orbit al angular momentu m / 
at every step of spin-up evolution, we used the lBeiger et al.l (|2010l ) re- 
sult: the Keplerian orbital frequency of a particle in the thin disk at a 
radius ro is well- approximated by the "Schwarzschildian/Newtonian" 
formula yGM/ro^ . This ansatz yields a surprisingly accurate deter- 
mination of I for ro — rms for a broad range of stellar masses and 
spin frequencies up to the mass-shedding limit. Here we show that 
this approximate approach allows for a reasonably accurate determi- 
nation of the rms as well. It corresponds to the d//dr = condition; 
the inner-edge of an accretion disk ro is a solution of 

where rcor is the corotation radius, Vm = (GM)~-'^/^M~2/7^4/7 jg 
the magnetospheric radius and ^ = ^m/^cor- Function /ms intro- 
duce d by Beiger et al. (201 J; reducing to /ms = 1 in the Eq. 17 of 
iKluzniak &: RaDDaDorti«2007l ) ensures a proper behavior of the torque 
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